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ABSTRACT 


Closed-form  expressions  for  the  perturbation  of  the  surface-wave  spectrum 
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ly  perturbed  spectral  component  as  a  function  of  these  parameters.  The  per¬ 
turbation  due  to  more  complicated  features  may  be  obtained  from  the  results 
presented  in  this  work  through  Fourier  expansion  techniques. 
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INTRODUCTION 


The  purpose  of  this  study  is  to  examine  quantitatively  the  perturbation  in 
the  surface-wave  height  spectrum  due  to  the  interaction  of  the  surface  waves 
with  the  current  field  of  an  internal  wave.  The  problem  has  been  studied  else¬ 
where  (e.g.,  Refs.  1  and  2  and  references  cited  therein),  but  we  feel  that  there 
are  enough  subtleties  to  warrant  further  examination.  In  particular,  we  want 
to  estimate  the  component  in  the  surface-wave  spectrum  that  is  most  strongly 
perturbed  for  a  given  internal-wave  phase  speed  and  wavelength. 


FORMULATION 


We  assume  that  our  current  field  is  given  by 

U  =  -U0  sin  [k,(x  -  C/t)l,  ^  x  -  C,t  < 

k,  k, 

=  0  elsewhere. 


(1) 


Equation  1  describes  a  sinusoidal  feature  one  wavelength  ( =  2 x/k,)  in  extent, 
propagating  along  the  x  axis  with  phase  velocity  C,.  This  simplified  represen¬ 
tation  of  the  internal-wave  current  field  contains  the  features  necessary  for  our 
purposes  and  keeps  the  mathematics  relatively  uncomplicated.  A  more  realis¬ 
tic  internal-wave  current  field  could  be  given  by  a  linear  superposition  of  terms 
such  as  that  of  Equation  1 . 

We  have  shown3  that  the  ratio  of  the  perturbed  surface-wave  height  spec¬ 
trum  S(k;  x,  t)  at  position  x  and  time  t  in  a  general  two-dimensional  current 
field  to  the  equilibrium  spectrum  $„( k)  is  given  by 


S(k;x,f) 


(2) 


where 


S,„(k)  1  +  P(k;x,l)  ’ 

P(k;x,r)  =  N"( k)  kj  —p1  ]  exp  [-  0(k')rff']  df ,  (3) 


and  k)  is  the  equilibrium  action  spectrum.  In  Eq.  3,  the  time  dependence 
of  x  and  k  is  determined  by  the  coupled  equations, 


dx, 

~dt 


dwo 

Ik, 


+  U, 


(4a) 


'B.  A.  Hughes,  Speckle  Noise  and  Pattern  Detection  in  SAR  Imagery  of  Internal  Waves,  81-12, 
Defense  Research  Establishment  of  the  Pacific  (Dec  1981). 

2D.  R.  Thompson  and  R.  F.  Gasparovic,  "Intensity  Modulation  in  SAR  Images  of  Internal  Waves” 
(to  be  published  in  Nature,  Apr  1986). 

5D.  R.  Thompson,  "Intensity  Modulation  in  Synthetic  Aperture  Images  of  the  Ocean  Surface  and 
the  Wave-Current  Interaction  Process,"  Johns  Hopkins  APL  Tech.  Dig.  6,  J46  (1985). 


1 


THE  JOHNS  HOPKINS  UNIVERSITY 

APPLIED  PHYSICS  LABORATORY 

LAUREL.  MARYLANO 


and 


dk,  ,  dUj 

~dt  ~  ~ j  ax7 


(4b) 


where  U,  is  the  rth  component  of  the  current  Held  and  «0  is  the  intrinsic 
surface-wave  frequency  given  by  the  gravity-capillary  dispersion  relation, 


wo  =  fM 


Also  in  Eqs.  3  and  4,  =  1  /N^  (k)  and  repeated  indices  are  summed. 

The  function  /3(h)  in  Eq.  3  is  the  so-called  relaxation  rate,  and  it  depends  on 
wind  velocity  as  well  as  k.  Later,  we  shall  see  that  the  magnitude  of  the  spec¬ 
tral  perturbation  is  quite  sensitive  to  /3(k). 

The  equations  for  the  spectral  perturbation  given  above  are  general,  and  their 
solution  can  be  rather  complicated.  We  now  want  to  simplify  our  particular 
problem  further  by  assuming  that  the  internal-wave  currents  are  always  much 
less  than  the  group  velocity  (  =  du)0/dk,)  of  the  surface  waves.  The  group  ve¬ 
locity  for  gravity-capillary  waves  has  a  minimum  value  of  about  18  cm/sec 
at  a  wavelength  of  around  4  cm.  Therefore,  if  we  require  U0  in  Eq.  1  to  be 
about  1  cm/sec  or  less,  we  can  neglect  £/,  in  Eq.  4a.  If  we  further  require  the 
current  gradient  to  be  small  compared  to  l//3(k),  we  see  from  Eq.  4b  that  the 
percent  change  in  k  over  the  relaxation  time  (=  1//3)  will  also  be  small.  For 
1-m  waves,  this  requires  the  current  gradient  (  =  k,U0)  to  be  less  than  around 
10'3  sec"1.  With  UQ  *  1  cm/sec,  this  requires  that  k,  s  0.1  m"1  or  the  wave¬ 
length  of  the  internal-wave  feature  to  be  greater  than  about  50  m. 

With  the  internal-wave  current  given  by  Eq.  1  and  the  restrictions  on  U0 
and  k,  discussed  above,  Eq.  3  reduces  to 


fllnAL 


Villi  Y  j.  .j  h 

PQnx,t)  =  k,U0kx  — — 5  exp(-/3f)  exp(/3 t')  cos  [  MC„  ~  C,)t'\  dt'  , 

OKx  J  rft 

°  (6) 


where  kx  and  denote,  respectively,  the  components  of  the  surface-wave 
number  and  group  velocity  along  the  x  axis.  In  writing  Eq.  6,  we  have  used  the 
relationship  x=  C^t  obtained  by  integrating  Eq.  4a  and  neglecting  the  (/,  term. 
In  a  coordinate  system  fixed  to  the  moving  current  feature,  this  relationship 
becomes 


I 


I 


x  =  (C^  -  C,)t  , 


(7) 


where  we  have  assumed  that  the  origin  of  the  two  systems  coincides  at  t 
We  may  rewrite  Eq.  6  in  terms  of  x  using  Eq.  7  to  obtain 


fl(k/) 


k,U0  din (  0x  \ 
c^-c,  x  dkx  expV~qK-  C,) 


mirtfjr.T/*/ ) 
-w/kf 


cos  k,x'  dx' 


=  0. 


(8a) 


2 
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for  Cg  z  C,  and  x  z  -  x/k,,  and 

k,U0  dlnN„  (  0x  \ 

'M  =  c^Q  **  inf  a>(-c^c,) 


JmaxCx, -»/*/)  /  qx,  \ 

exp(  )  cos  k,x'  dx' 

t/kj  ' 


(8b) 


for  Cg  <  C/  and  x  <,  x/k,.  The  integration  limits  span  the  region  of  nonzero 
current  and  depend  on  the  sign  of  Cg  -  C,,  which  determines  whether  x  in¬ 
creases  or  decreases  with  time.  Equations  8  may  be  integrated  to  obtain 


31nN„ 

P(k^c)  =  k,U0  k.  * 


1 


9kx  [0 2  +  k1,  (C„  -  Ci)1] 


2-1  ‘A 


^cos  (k,x  -  5)  +  cos  S  expj^-  ^/x  +  s%n  *)  j  j  ^  (9a) 


with 


and 


tan  6 


*;(C„  -  C,) 

0 


sgn 


1  for  Cgg  z  C, 
-1  for  Cfx  <  C, 


(9b) 


(9c) 


for  -x/k,  <  x  <  x/k, .  For  x  >  x/k, ,  we  Find 


/'(k^r)  =  -  (1  +  sgn)  k,U0  kx 


din Ngf 

~JkT 


0 

02  +  tlCg,  -  C,)1 


•  sinh 


(9d) 


for  x  <  -x/k,,  P(k^ir)  equals  the  result  in  Eq.  9d  with  the  factor  -  (1  +  sgn) 
replaced  by  (1  -  sgn).  Finally,  for  x  <  -x/k,  and  Cg  z  C,  or  for  x  a  x/k, 
and  Cg  <  C/,  PQl,  x)  =  0. 
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DISCUSSION 


Equations  9  are  the  desired  expressions  for  the  perturbation  of  the  kth  com¬ 
ponent  of  the  surface-wave  spectrum  as  a  function  of  position  along  the  current 
feature.  Note  that  if  =  C,,  these  equations  reduce  to 

31nN„  1 

AM)  =  k,U0  kx  —  cos  k,x  (10a) 

dkx  p(k) 


for  |x|  <  r/k\  and  zero  otherwise.  For  a  Pierson-Moskowitz  height  spectrum 
with  a  cos"  <t>/2  angular  dependence  where  <t>  is  the  angle  between  k  and  the  wind 
direction,  we  find  that 


din Af„ 

*  dkx 


cos  0  + 


n  <t> 

-tan  -  sin  0 
2  2 


cos  0  , 


(10b) 


where  g  is  the  acceleration  of  gravity,  Vw  is  the  wind  speed,  and  0  is  the  angle 
between  k  and  the  x  axis  (internal-wave  propagation  direction).  Inserting  this  ex¬ 
pression  into  Eq.  10a,  we  obtain 


AM) 


179 

1.48g2N 

„  n  *  ■  J 

1  dU 

L\2 

Kk1) 

cos  0  —  tan  -  sin  0  cos  0 

2  2  J 

0(k)  dx 

(10c) 


where  dU/dx  is  the  current  gradient.  For  wind  speeds  of  about  3  m/sec  or  more 
and  k  values  corresponding  tc  waves  of  1  m  or  less  propagating  in  the  wind  direc¬ 
tion  (<t>  =»  0°),  Eq.  10c  can  be  reduced  further  to 

9,1  dU 

AM)  =  -  cos2  0  -  —  ,  (lOd) 

2  0(k)  dx 

an  equation  that  has  been  used  to  estimate  the  modulation  observed  in  Seasat 
synthetic  aperture  radar  (SAR)  images. 

From  Eqs.  10,  we  see  that  the  perturbation  varies  directly  as  the  current  gra¬ 
dient  and  inversely  with  the  relaxation  rate  /3(k).  However,  it  is  important  to 
remember  that  Eqs.  10  hold  only  for  the  so-called  resonance  condition  when  the 
x  component  of  the  group  velocity  for  surface  waves  of  wave  number  k  is  equal 
to  the  phase  speed  of  the  internal-wave  feature.  Thus,  for  a  given  k  value  (e.g., 
determined  by  the  frequency  of  a  radar  transmitter),  this  resonance  condition 
will  not,  in  general,  be  met,  and  the  expressions  given  by  Eqs.  9  must  be  used 
to  compute  the  spectral  perturbation.  In  fact,  we  will  show  in  the  following  dis¬ 
cussion  that  the  maximum  perturbation  does  not  usually  occur  for  the  k  value 
that  yields  the  resonance  condition. 

In  order  to  understand  how  the  perturbation  given  by  Eqs.  9  varies  with  k 
for  a  given  current  feature,  it  is  useful  to  first  see  how  0(k)  and  Cg  vary  with 
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0  40  80  120  160  200 

k  (rad/m) 


Figure  1  /(-dependence  of  the  surface-wave  relaxation 
rate,  0,  and  group  velocity.  The  0  curves  are  for  wind 
speeds  V„  of  3  and  6  m/sec  and  for  wind  direction 
with  respect  to  k,  <t>  =  0*  and  <t>  =  60°.  The  parameteri¬ 
zation  of  /3<k)  is  that  of  Hughes.4 


k.  As  discussed  in  Refs.  3  and  4,  the  behavior  of  /3(k)  is  a  topic  of  current  re¬ 
search.  However,  we  have  found  recently  that  the  full  solution  of  Eqs.  3  and 
4  with  /3(k)  given  by  Hughes4  yields  good  agreement  with  measured  spectral  per¬ 
turbation  data  from  the  Georgia  Strait  Experiment.5  Therefore,  we  show  the  k 
dependence  of  Hughes’  0(k)  in  the  upper  plot  of  Fig.  1  for  wind  speeds  of  3 
and  6  m/sec  (measured  10  m  above  the  surface)  blowing  along  the  direction  of 
k  and  for  3  m/sec  blowing  at  an  angle  <t>  =  60°  from  the  k  direction.  Note  that 
for  all  cases  shown  /3(k)  increases  slowly  for  k  a  60  m'1  but  decreases  rather 
rapidly  for  A:  s  60  m'1.  The  surface-wave  group  velocity,  shown  in  the  lower 
plot  of  Fig.  1 ,  is  just  the  derivative  with  respect  to  k  of  the  gravity-capillary  dis¬ 
persion  relation  given  by  Eq.  5.  It  is,  in  fact,  almost  constant  at  about  0.2  m/sec 
between  k  =  60  and  200  m"1  but  varies  like  k~'A  for  k  ■£>  60  m'1  (the  gravity- 
wave  region). 

Thus  we  see  that  for  60  <  k  <  200  m'1,  the  entire  k  dependence  in  P(k^c) 
is  essentially  due  to  that  of  0(k).  Furthermore,  for  internal-wave  wavelengths  great¬ 
er  than  about  60  m  (k,  s  0.10  m'1)  and  -  0.1  <  C,<,  0.5  m/sec,  k,  ■  ( Cg  - 
C,)  is  small  compared  to  0  for  roughly  30  <  k  <  200  m~'.  For  these  condi¬ 
tions,  we  have  essentially  the  resonance  condition,  and  Eqs.  10  for  PQiyX)  ap¬ 
proximate  closely  the  more  general  expressions  of  Eqs.  9. 

On  the  other  hand  for  k  <,  30  m  '1,  the  k  dependence  of  both  /3(k)  and  Cg 
must  be  considered  in  estimating  the  maximum  perturbation.  In  Fig.  2,  we  show 


4B.  A.  Hughes,  “The  Effect  of  Internal  Waves  on  Surface  Wind  Waves;  Theoretical  Analysis," 
J.  Geophys.  Res.  83C,  455  (1978). 

5D.  R.  Lyzenga  and  R.  A.  Schuehman,  eds..  Interim  Report  on  the  Georgia  Strait  Experiment, 
Vol.  I,  Environmental  Res.  Inst,  of  Michigan  (Dec  1984). 
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Figure  2  Amplitude  factor  as  a  function  of  k  for  various  internal-wave  phase 
speeds.  The  four  curves  (solid,  dash,  short  dash,  and  dashfdot)  in  each  of  the 
plots  are  computed  for  internal-wave  wavelengths  of  60,  80,  100,  and  120  m, 
respectively.  The  k  vector,  internal-wave  propagation  direction,  and  wind  direc¬ 
tion  are  colinears  (6  =  0°,  <t>  =  0°),  and  the  wind  speed  is  3  m/sec  for  all  cases. 


plots  of  the  amplitude  factor  in  Eq.  9a  (the  term  that  multiplies  the  |  |  brackets) 
as  a  function  of  k  for  Fixed  internal-wave  phase  speeds  C, .  The  results  for  var¬ 
ious  wave  numbers  k,  corresponding  to  internal-wave  wavelengths  of  60,  80, 
100,  and  120  m  are  shown  by  the  solid,  dash,  short  dash,  and  dash/dot  curves, 
respectively.  The  arrows  on  the  k  axis  of  the  plots  in  Fig.  2  mark  the  k  value 
for  whic'n  the  resonance  condition  (Cgx  =  C,)  holds.  All  of  the  curves  in  the 
figure  were  computed  using  the  /3(k)  corresponding  to  a  3  m/sec  wind  along  the 
direction  of  k  (<£  =  0°),  and  with  the  angle  between  k  and  the  internal-wave 
propagation  direction  equal  to  zero  (i.e.,  6  =  0  and  Cgx  =  CK).  Also,  we  have 
assumed  for  these  and  all  subsequent  computations  that  the  peak  current  gra¬ 
dient,  krU0  =  101  sec'1,  remains  constant  and  that  the  kxd\rrfJeq/dkx  term  in 
Eqs.  9  is  given  by  Eq.  10b  with  n  =  4. 

We  can  see  from  Fig.  2  that  for  30  <  k  <  200  m  1  the  amplitude  factor  is 
nearly  independent  of  the  properties  of  the  internal  wave  (for  the  range  of  k, 
and  C,  values  considered)  as  expected  from  our  previous  discussion,  and  its  vari¬ 
ation  over  this  range  is  due  entirely  to  the  behavior  of  /3(k).  For  k  s  30  m  ', 
we  begin  to  see  the  dependence  of  the  amplitude  factor  on  kt  and  C, .  In  par¬ 
ticular,  we  note  that  for  a  given  C,,  the  largest  amplitude  occurs  for  the  smallest 
k,  and  that  the  largest  amplitudes  overall  occur  for  the  largest  phase  speed 
C,  =0.5  m/sec.  This  behavior  is  simply  a  result  of  the  k  dependence  of  0  and 
Cg  in  the  denominator  of  the  amplitude  term  in  Eq.  9a.  Since  0(k)  increases  with 
k  faster  than  Ct  decreases,  this  denominator  will  always  have  a  minimum  at  a 
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Figure  3  Amplitude  factor  as  a  function  of  k  for  various  angles  6  between  the 
k  direction  and  the  internal-wave  propagation  direction.  The  four  curves  (solid, 
dash,  short  dash,  and  dash/dot)  in  each  of  the  plots  are  computed  for  internal- 
wave  wavelengths  of  60, 80, 100,  and  120  m,  respectively.  The  internal-wave  phase 
speed  is  0.2  m/sec,  and  the  wind  speed  is  3  m/sec  along  the  k  direction  (<f>  = 
0°)  for  all  cases. 


k  value  less  than  that  for  which  Cg  =  C,.  The  value  of  this  minimum  (which 
produces  a  maximum  in  P(k ,x))  is,  of  course,  controlled  by  the  value  of  k,  and 
C,  for  the  particular  internal  wave  in  question.  A  decrease  in  k,  or  an  increase 
in  C,  will  lower  the  value  of  this  minimum  and  shift  its  position  to  a  smaller 
k  value. 

In  Fig.  3,  we  show  how  the  amplitude  factor  varies  as  a  function  of  6,  the 
angle  between  k  and  the  internal-wave  propagation  direction.  In  each  of  the  plots 
in  this  Figure,  we  have  assumed  an  internal-wave  phase  speed  of  C,  =  0.2  m/sec 
and  a  (3(k)  corresponding  to  3  m/sec  winds  along  the  direction  of  k  (i.e.,  <t>  = 
0°),  as  before.  The  solid,  dash,  short  dash,  and  dash/dot  curves  in  each  plot 
are  again  for  internal-wave  wavelengths  of  60,  80,  100,  and  120  m,  respectively, 
and  the  arrows  indicate  the  k  value  where  the  resonance  condition  holds.  (There 
can  be  no  resonance  condition  for  0  >  90°. )  One  can  see  from  the  plots  in  this 
figure  (and  the  C,  =  0.2  m/sec  plot  in  Fig.  2,  which  is  for  0  =  0°)  that  the 
amplitude  factor  decreases  as  0  approaches  90°  due  to  the  cos2  0  dependence 
of  kxd\nNeq/dkx  (as  given  by  Eq.  10b  with  <t>  =  0°).  As  0  becomes  larger  than 
90°,  the  k 2  (C^  -  C,)1  term  in  the  denominator  of  Eq.  9a  begins  to  grow.  It 
can  be  seen  from  Fig.  3  that  the  size  of  this  term  along  with  the  cos2  0  depen¬ 
dence  controls  the  behavior  of  the  amplitude  factor  for  S  i  120°. 

Figure  4  illustrates  how  the  wind-direction  dependence  of  0(k)  can  affect  the 
amplitude  factor.  In  each  of  the  plots  in  this  figure,  we  have  again  chosen  an 
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Figure  4  Amplitude  factor  as  a  function  of  k  for  various  angles  between  the 
k  direction  and  wind  direction.  The  four  curves  (solid,  dash,  short  dash,  and 
dash/dot)  in  each  of  the  plots  are  computed  for  internal-wave  wavelengths  of 
60, 80, 100,  and  120  m,  respectively.  The  internal-wave  phase  speed  is  0.2  m/sec, 
the  wind  speed  is  3  m/sec,  and  the  k  direction  is  along  the  internal-wave  propa¬ 
gation  direction  (0  =  0*)  for  all  cases. 


internal-wave  phase  speed  of  Q  =  0.2  m/sec  and  a  wind  speed  of  3  m/sec  as 
in  Fig.  3.  For  this  case,  however,  we  have  assumed  that  the  k  vector  is  pointing 
along  the  x  axis  ( d  =  O’)  for  each  plot  and  that  the  wind  direction  is  at  an  angle 
<t>  with  respect  to  k  as  indicated.  As  before,  solid,  dash,  short  dash,  and  dash/dot 
curves  in  each  plot  represent  internal-wave  wavelengths  from  60,  80,  100,  and 
120  m,  respectively.  The  resonance  condition  occurs  at  k  =  75  m"1.  Since  /3(k) 
decreases  like  cos  <t>  (see  Fig.  1),  we  see  from  Fig.  4  that  the  amplitude  factors 
increase  with  increasing  #  as  expected.  In  particular,  it  can  be  seen  that  the  k 
value  corresponding  to  the  maximum  amplitude  factor  increases  with  increasing 
tt>  (decreasing  /3).  For  example,  as  <p  increases  from  0  to  60°,  a  comparison  of 
the  internal-wave  wavelength  equals  100  m  curve  in  the  C,  =  0.2  m/sec  plot 
in  Fig.  2  (for  which  <t>  =  0°)  with  the  corresponding  curve  in  the  <t>  =  60°  plot 
in  Fig.  4  shows  that  the  position  of  maximum  amplitude  moves  from  10  to  about 
18  m'1  and  increases  in  magnitude  from  0.2  to  0.3. 

In  order  to  see  how  the  perturbation  varies  over  the  extent  of  the  internal- 
wave  feature,  we  show  in  Fig.  5  plots  of  the  assumed  internal-wave  current  (giv¬ 
en  by  Eq.  1)  and  the  resulting  spectral  perturbation  P( k,  x)  (given  by  Eqs.  9) 
for  specified  k  values  as  a  function  of  x.  In  constructing  the  plots  shown  in  Fig. 
5,  we  have  chosen  an  internal-wave  wavelength  of  100  m,  a  phase  speed  of  0.2 
m/sec,  and  a  wind  speed  of  3  m/sec  along  the  propagation  direction.  The  wave 
number  vector  of  the  surface  waves  is  chosen  to  be  30  and  10  m'1  in  the  mid- 
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Figure  5  Current  and  corresponding  perturbation  as 
a  function  of  position  for  k  =  30  m'1  (middle  plot)  and 
10  m*1  (lower  plot).  The  Internal-wave  phase  speed  is 
0.2  m/sec;  the  wind  speed  is  3  m/sec;  and  the  k  vec¬ 
tor,  internal-wave  propagation  direction,  and  wind 
direction  are  colinear  (6  =  O',  4  =  O')  for  each  case. 

The  •  symbol  shows  the  perturbation  predicted  by  the 
resonance  approximation. 

die  and  lower  plots,  respectively,  with  both  vectors  pointing  in  the  internal-wave 
propagation  direction.  (These  conditions  correspond  to  0  =  0°  and  =  0°.) 
The  •  symbols  in  the  lower  two  plots  show  the  perturbations  predicted  by  the 
resonance  approximation  given  by  Eq.  10c.  We  see  from  Fig.  5  that,  for  k  = 
30  m*1,  the  variation  of  the  perturbation  as  a  function  of  position  in  the  internal- 
wave  current  feature  given  by  the  resonance  condition  agrees  well  with  that  giv¬ 
en  by  the  more  general  expressions  of  Eqs.  9.  In  particular,  the  perturbation  for 
this  case  is  roughly  proportional  to  the  gradient  of  the  current  except  near  the 
current  boundaries  where  the  transient  terms  in  Eqs.  9  are  important.  On  the 
other  hand,  for  the  k  =  10  m"1  case  shown  in  the  bottom  plot  of  Fig.  5,  we 
see  that  the  maximum  perturbation  is  increased  by  about  a  factor  of  two  over 
the  k  =  30  m'1  case.  This  increase  is  smaller  by  a  factor  of  about  two  than  the 
increase  predicted  by  the  resonance  approximation  (and  shown  by  the  •  symbols 
in  the  plot),  and  the  position  of  the  maximum  perturbation  is  at  x  *  15  m  rath¬ 
er  than  x  =  0  m  due  to  the  phase  shift  6  that  appears  in  Eqs.  9.  Note  also  that 
for  the  it  =  10  m'1  case,  there  is  a  non-negligible  perturbation  extending  for 
about  20  m  or  so  beyond  the  current  boundary  at  x  =  50  m,  a  feature  also  not 
predicted  by  the  resonance-condition  computations. 
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Figure  6  Current  and  corresponding  perturbation  as 
a  function  of  position  for*  =  30  m*1  (middle  plot)  and 
10  m"1  (lower  plot).  The  internal-wave  phase  speed  is 
0.2  m/sec;  the  wind  speed  is  6  m/sec;  and  the  k  vec¬ 
tor,  internal-wave  propagation  direction,  and  wind 
direction  are  collnear  ( 6  =  0 ",</>  =  0*)  for  each  case. 
The  •  symbol  shows  the  perturbation  predicted  by  the 
resonance  approximation. 


Figure  6  was  computed  in  exactly  the  same  manner  as  Fig.  5  except  that  a 
wind  speed  of  6  m/sec  (again  in  the  internal-wave  propagation  direction)  was 
assumed  in  the  computation  of  0(k).  The  important  feature  here  is  that  the  reso¬ 
nance  approximation  gives  results  that  are  in  closer  agreement  with  the  more  gener¬ 
al  computation  than  were  found  for  the  3  m/sec  wind  speed.  One  can  see  that, 
for  the  6  m/sec  wind,  the  relative  increase  in  the  maximum  perturbation  between 
k  =  10  and  30  m'1  is  about  a  factor  of  three.  Also  by  comparing  Figs.  5  and 
6,  we  see  that  the  absolute  magnitude  of  the  perturbation  (computed  using  Eqs. 
9)  decreases  by  only  a  factor  of  about  1 .3  in  going  from  3  to  6  m/sec  winds. 
This  is  a  smaller  reduction  than  the  factor  of  two  predicted  by  the  resonance 
approximation. 

We  illustrate  the  difference  in  P(k,  x)  for  0  =  0"  and  9  =  180*  in  Fig.  7, 
where  the  middle  and  lower  plots  show  the  comparison  for  k  =  30  m~l  and  k 
=  10  m'1,  respectively.  An  internal-wave  phase  speed  of  C,  =  0.2  m/sec  and 
a  wavelength  of  100  m  were  again  assumed  for  these  computations,  along  with 
a  wind  speed  of  3  m/sec  in  the  k  direction  (</>  =  0").  The  k  =  30  m'1  results 
are  rather  similar  both  in  shape  and  in  magnitude.  Those  for  k  =  10  m'1  are 
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Figure  7  Current  and  corresponding  perturbation  as 
a  function  of  position  for  k  =  30  m' 1  (middle  plot)  and 
10  m'1  (lower  plot).  The  internal-wave  phase  speed  is 
0.2  m/sec  and  the  wind  speed  is  3  m/sec  along  the  fr 
direction^  =  0°)  for  each  case.  The  solid  curve  in  each 
of  the  perturbation  plots  is  for  Ar  along  the  internal-wave 
propagation  direction  ( 6  =  0°),  while  the  dashed  curve 
is  for  k  opposite  the  propagation  direction  (6  =  180*). 
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roughly  130°  out  of  phase,  with  the  maximum  at  x  *  15  m  for  the  8  =  0°  curve 
being  about  a  factor  of  two  larger  than  the  maximum  at  x  *  -  20  m  for  the 
0  =  180°  curve.  This  figure  illustrates  the  fact  that  the  0  dependence  of  fl(k, 
x)  is  more  complicated,  especially  at  lower  k  values,  than  the  simple  cos2  0  de¬ 
pendence  given  by  the  resonance  approximation  of  Eq.  10c. 

As  a  practical  application  for  the  results  presented  above,  we  point  out  that 
in  the  context  of  Bragg  scattering  theory  the  intensity  modulation  observed  in 
the  radar  image  of  a  particular  current  feature  is  expected  to  be  proportional 
to  the  perturbation  in  the  surface- wave  spectrum  at  the  Bragg  wavelength.  (For 
a  discussion  of  Bragg  scattering  of  electromagnetic  radiation  from  a  rough  sur¬ 
face,  see  Ref.  6.)  For  an  L-band  SAR,  for  example,  the  Bragg  wavelength  is  in 
the  0.2  m  range,  which  corresponds  to  Ar  *  30  m"1.  Therefore,  if  such  a  radar, 
looking  along  the  propagation  direction,  were  to  image  our  nominal  internal- 
wave  current  feature  of  wavelength  100  m  and  phase  speed  0.2  m/sec,  and  if 


6J.  W.  Wright,  "Backscattering  from  Capillary  Waves  with  Application  to  Sea  Clutter,”  IEEE 
Trans.  Antennas  Propag.  AP-14,  749  (1966). 


11 


v_ 


THE  JOHNS  HOPKINS  UNIVERSITY 

APPLIED  PHYSICS  LABORATORY 

LAUREL.  MARYLAND 


-100  -75  -50  -25  0  25  50  75  100 


x  (m) 

Figure  8  Current  and  corresponding  perturbation  as 
a  function  of  position  for  a  nominal  k  =  30  m'1  and 
k  optimized  for  maximum  perturbation  for$  =  0*  (mid¬ 
dle  plot)  and  <t>  =  60*  (lower  plot).  The  internal-wave 
phase  speed  Is  0.2  m/sec,  the  wind  speed  is  3  m/sec, 
and  the  k  direction  is  along  the  propagation  direction 
(6  =  0“)  for  each  case. 


the  wind  speed  were  3  m/sec  also  along  the  propagation  direction,  we  would 
expect  to  see  an  image  intensity  modulation  across  the  current  feature  that  varies 
like  the  solid  curve  in  the  middle  plot  of  Fig.  8.  However,  from  our  previous 
discussion  we  know  that  the  degree  of  modulation  for  an  internal-wave  feature 
of  this  type  is  substantially  enhanced  at  lower  k  values.  In  particular,  we  see  from 
the  phase  speed  =  0.2  m/sec  plot  in  Fig.  2  that  the  maximum  amplitude  factor 
for  a  100  m  feature  occurs  at  about  k  =  10  m  l.  Thus,  by  “tuning”  our  radar 
to  a  wavelength  of  about  0.63  m  (which  corresponds  to  k  =  10  m  '),  we  would 
expect  to  see  the  intensity  modulation  given  by  the  dashed  curve  in  the  middle 
plot  of  Fig.  8.  By  tuning  to  this  “optimum  frequency,”  we  obtain  a  factor  of 
two  gain  in  the  maximum  perturbation  and  also  increased  modulation  (brighter, 
since  F(k,  x)  <  0)  for  x  >  50  m  outside  the  current  feature  as  shown  in  this  plot. 

The  bottom  plot  in  Fig.  8  shows  a  similar  example  for  which  the  wind  is  as¬ 
sumed  to  be  blowing  at  60°  to  the  propagation  direction  and  radar  look  direc¬ 
tion  (i.e.,  8  =  0'  and  <t>  =  60’).  The  solid  curve  in  this  plot  shows  the  expected 
intensity  modulation  for  =  30  m1,  while  the  dashed  curve  shows  the  modu¬ 
lation  for  the  optimum  k  =  18  m  "1  obtained  from  the  <t>  =  60°  plot  in  Fig.  4. 
One  can  see  that  the  k  =  30  m'1  modulation  is  already  larger  than  the  cor- 
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responding  =  0*  case  due  to  the  decrease  in  0(k)  with  increasing  However, 
an  additional  factor  of  1 .5  can  be  gained  in  the  maximum  modulation  by  tuning 
to  k  =  18  m'1.  Also,  as  was  seen  for  the  4>  =  0°  case,  this  tuning  yields  a  sub¬ 
stantial  increased  brightening  outside  the  current  feature  for  x  >  SO  m. 

We  conclude  that  for  given  wind  conditions,  internal-wave  properties,  and  view¬ 
ing  geometry,  one  can,  at  least  in  principle,  tune  a  radar  as  illustrated  in  the  above 
examples  to  obtain  the  optimum  intensity  modulation  in  the  image.  Conversely, 
it  may  also  be  possible  to  use  a  knowledge  of  the  return  signal  (along  with  view¬ 
ing  geometry  and  wind  conditions)  to  deduce  the  properties  of  the  current  feature. 


SUMMARY  AND  CONCLUSIONS 


In  this  study,  we  have  developed  closed-form  expressions  for  the  perturbation 
of  a  two-dimensional  surface-wave  height  spectrum  due  to  a  one-dimensional 
sinusoidal  current  feature  with  specified  wavelength  and  phase  speed.  These  ex¬ 
pressions  were  derived  from  a  general  formulation  of  the  wave-current  interac¬ 
tion  problem,3  with  the  restriction  that  only  currents  whose  magnitude  is  less 
than  a  few  centimeters  per  second  and  whose  gradient  magnitude  is  less  than 
a  few  10'3  sec-1  be  considered.  The  expressions  have  been  checked  against  the 
general  calculation  for  a  number  of  different  cases  (which  satisfy  the  above  re¬ 
strictions),  and  excellent  agreement  was  found  in  each  case. 

We  have  used  these  expressions  to  examine  how  the  magnitude  of  the  spectral 
perturbation  depends  on  quantities  such  as  the  surface-wave  wave  number  and 
propagation  direction,  internal-wave  phase  speed  and  wavelength,  and  wind  speed 
and  direction.  The  study  showed  that  regimes  exist  where  the  perturbation  is  sen¬ 
sitive  to  changes  in  each  of  these  parameters.  We  have  also  found  that  the  so- 
called  resonance  approximation  for  the  spectral  perturbation,  which  is  exact  when 
the  surface-wave  group  velocity  is  equal  to  the  phase  velocity  of  the  current  fea¬ 
ture,  is  reasonably  valid  for  surface-wave  numbers  greater  than  about  30  m'1 
with  wind  speeds  greater  than  about  3  m/sec.  However,  as  the  magnitude  of  the 
wave  number  decreases  below  30  m'1,  the  resonance  approximation  breaks  down 
rather  quickly,  and  the  more  general  expressions  must  be  used  to  obtain  accurate 
estimates  of  the  perturbation. 

The  behavior  of  the  perturbation  as  a  function  of  position  across  the  current 
feature  has  also  been  examined  as  a  function  of  the  various  parameters,  and  it 
was  found  that  the  position  of  the  maximum  and  minimum  values,  for  example, 
can  be  quite  sensitive  to  changes  in  these  parameters.  In  particular,  we  have  seen 
that,  for  the  combination  of  parameters  that  yield  the  largest  magnitude  pertur¬ 
bation,  the  position  of  the  maximum  (or  minimum)  is  usually  not  at  the  position 
of  maximum  (or  minimum)  current  gradient.  This  is  in  contrast  to  the  resonance 
approximation  that  predicts  that  the  perturbation  should  vary  like  the  current 
gradient.  Furthermore,  it  was  found  that,  in  many  cases,  significant  perturba¬ 
tion  can  exist  well  beyond  the  extent  of  the  current  feature— again  in  contradic¬ 
tion  to  the  resonance  approximation. 

Finally,  we  have  shown  how  it  may  be  possible,  at  least  in  the  context  of  Bragg 
scattering  theory,  to  tune  the  frequency  of  a  radar  using  our  results  in  order  to 
maximize  the  observed  intensity  modulation  from  a  given  current  feature.  On 


13 


REFERENCES 


1 B.  A.  Hughes,  Speckle  Noise  and  Pattern  Detection  in  SAR  Imagery  of  Inter¬ 
nal  Waves,  81-12,  Defense  Research  Establishment  of  the  Pacific  (Dec  1981). 

2D.  R.  Thompson  and  R.  F.  Gasparovic,  “Intensity  Modulation  in  SAR  Im¬ 
ages  of  Internal  Waves”  (to  be  published  in  Nature,  Apr  1986). 

3D.  R.  Thompson,  “Intensity  Modulation  in  Synthetic  Aperture  Images  of  the 
Ocean  Surface  and  the  Wave-Current  Interaction  Process,”  Johns  Hopkins  APL 
Tech.  Dig.  6,  346  (1985). 

4B.  A.  Hughes,  “The  Effect  of  Internal  Waves  on  Surface  Wind  Waves;  The¬ 
oretical  Analysis,”  J.  Geophys.  Res.  83C,  455  (1978). 

SD.  R.  Lyzenga  and  R.  A.  Schuehman,  eds.,  Interim  Report  on  the  Georgia 
Strait  Experiment,  Vol.  I,  Environmental  Res.  Inst,  of  Michigan  (Dec  1984). 

®J.  W.  Wright,  “Backscattering  from  Capillary  Waves  with  Application  to  Sea 
Clutter,”  IEEE  Trans.  Antennas  Propag.  AP-14,  749  (1966). 


V. 

9 


14 


>>>  kvc  ‘/vsi  w.  -rm  vv.*  mi-  o. 


INITIAL  DISTRIBUTION  EXTERNAL  TO  THE  APPLIED  PHYSICS  LABORATORY 

The  work  reported  in  TG  1356  was  done  under  Navy  Contract  N00024-85*C-5301.  This  work  is  related  to  Task  ZDAO.  which  is  supported  by  Defense  Ad¬ 
vanced  Research  Projects  Agency  (STO). 


Requests  for  copies  of  this  report  from  DoD  activities  should  bo  directed  to  DTICfFDRB, 

Cameron  Station,  Bldg.  5,  Alexandria,  Virginia  22304-6145  using  OTIC  Form  1  and,  If  necessary,  OTIC  Form  55. 


